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Abstract 

It is demonstrated how to obtain the least biased description of the single 
particle spectra measured in all multiparticle production processes by using 
information theory approach (known also as MaxEnt approach). The case of 
e + e~ annihilation in hadrons process is discussed in more detail as an example. 
Comparison between MaxEnt approach and simple dynamical model based on 



Information theory is nowdays widely used in all branches of sciences to provide 
the least biased, most probable description of data for which we have only limited 
amount of information [3 . The information is defined by means of information 
entropy, which, in turn, can assume extensive (Shannon) or nonextensive (Tsallis) 
form, depending on circumstances (see IH1 El for details and further references) . 
In our case where we are interested in single particle rapidity distributions, dN/dy, 
maximization of such entropy (known as MaxEnt procedure [3 El 13) with constraint 
given by the energy conservation leads to the following form of the probability distri- 



bution to find one of the N produced particles of transverse mass /it = J ^ 2 + (pt) 2 
in the longitudinal phase space described by rapidity y such that energy of this 
particle is E = fiT coshy and its longitudinal momentum is pi = /iTsinhy: 



P(y) = = J-exp g (-/V// r coshy). 



This form is identical with that used in statistical models but now Z q and j3 q are no 
longer free parameters to be fitted when comparing with experimental data but in- 
stead are given by the normalization condition and energy conservation constraint 1 , 



dyp(y) = 1 and / dy fi T ■ cosh y ■ [p(y)] c - 



N 



(where ±Y m are maximal rapidities available in rest frame of hadronizing source, 
see |H1 IH fo r details). As is demonstrated in |U |S] with q > 1 (responsible for 
dynamical fluctuations, see [HI EH] and g-inelasticity K q (connected to true inelasticity 
K q defining a fraction of the original available energy W used for the production of 
secondaries, K q = n q /(3 — 2q), cf. jl]) as the only parameters we can describe all 
data for pp and pp collisions as well as recent data for nuclear collisions. Here we 
show in Fig. 1 that we can also fit fairly well e + e~ annihilation data for which, by 
definition, K q = 1 (the whole energy must be used for the production, there are no 
leading particles) and q remains the only free parameter. It turns out that in this 



case q < 1 2 . This fact indicates (see [TU]) that in e + e~ processes (where K q = 1) one 
cannot obtain equilibrium state because the corresponding temperature parameter 
depends now on energy: T — 1/ (3 q = Tq — (1 — q)E and system is highly influenced 
by the conservation of energy constraint. In Fig. 1 we have also shown attempts to 
fit the same data using the well defined sequential decay (cascading) hadronization 
developed in [7] instead of instantaneous hadronization 3 . Cascade here is defined by 
the sequence of decay processes, Mi -> [M^\ = kxMi] + [Afgi = k 2 Mi] (proceeding 
until Mi > 2/jt) in which all dynamics is summarily described by distribution P(k) 
of decay parameters fcj G (0, 1) (such that in each vertex k\ + k 2 < 1, see [7j for 
details). It is easy to realise that in this case dN/dy is given by a kind of random 
walk in rapidity space, which, in turn, results in its gaussian-like shape (cf. Fig. 1). 
However, one cannot find any reasonable values of decay parameters k\ to fit data 
with such distribution, whereas MaxEnt in its nonextensive version can do it with 
q = 0.6. The remaining small discrepancies for small and large rapidities indicate 
therefore a need for some additional information to be incorporated here (see, for 
example, discussion in 5J). When introduced it could then be again checked against 
experimental data. In this way we could always obtain good description of data 
with only minimal number of assumptions, i.e., with minimal information content 
[HUE]- This is the main advantage of the method presented here and in our opinion 
it deserves further detailed investigations by using wide spectrum of the available 

multiparticle production data of all kinds. 

2 To be contrasted with q > 1 values needed to describe the pr distributions in the same processes 
instead 0. 

3 Cascade could be regarded as a viable alternative because of the a priori cascade character of 
the quarks gluon and gluon — > quark-antiquark processes preceding hadronization expected to be 
present in e + e~ annihilations pQ. 
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Fig.l Upper-left panel: visualization of differences between sequential cascade hadronization (CAS) |7| and 
instantaneous one described by MaxEnt 0. Upper-right panel: closer look at the way dN/dy arises in CAS 
(resembling random walk in j/-space). Lower panels: results of calculations of dN/dy using CAS and MaxEnt and 
compared with e+e — annihilation data by ALEPH 0. Left panel - comparison with CAS model giving the same 
multiplicity |7| (it can be extremely well approximated by gaussian fit) and with MaxEnt for K q = 1 and q = 1 (it 
does not fit data). Right panel - gaussian fit has been replaced by MaxEnt fit with total inelasticity K q = 1 and 
with varying q; for q = 0.6 we obtain quite good agreement with data. 
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